Raman Photoassociation of Bose-Fermi Mixtures and the Subsequent Prospects for 

Atom-Molecule Cooper Pairing 
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We theoretically investigate Raman photoassociation of a degenerate Bose-Fermi mixture of atoms 
and the subsequent prospect for anomalous (Cooper) pairing between atoms and molecules. Stable 
fermionic molecules are created via free-bound-bound stimulated Raman adiabatic passage which, 
in contrast to purely bosonic systems, can occur in spite of collisions. With the leftover atomic 
condensate to enhance intrafermion interactions, the superfiuid transition to atom-molecule Cooper 
pairs occurs at a temperature that is roughly an order of magnitude below what is currently feasible. 
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Photoassociation occurs when two atoms absorb a 
laser photon thereby jumping from the free two- 
atom continuum to a bound molecular state. Neutral- 
atom statistics is determined by the number of neutrons 
in the nucleus-odd for fermions and even for bosons. 
Similarly, the sum of the total number of neutrons in 
the nuclei of the constituent atoms determines neutral- 
molecule statistics. Molecules formed by photoassocia- 
tion of two fermions will accordingly result in a boson, 
whereas fermionic molecules are born of a boson and a 
fermion. Given degenerate Bose-Fermi atoms, two ques- 
tions arise: Will the atoms photoassociate with into an 
arbitrary number of stable Fermi molecules? If so, is it 
possible to realize atom-molecule Cooper pairing? 

First introduced to explain superconductivity, anoma- 
lous quantum correlations between two degenerate elec- 
trons with equal and opposite momenta-Cooper pairs- 
are due physically to an electron-electron attraction me- 
diated by the exchange of lattice- vibration-generated 
phonons |2( , and are a manifestation of fermionic super- 
fluidity 3]. Anomalous pairing between different chem- 
ical species was immediately suggested to explain the 
larger excitation energy for nuclei with even rather than 
odd numbers of nucleons p|, although it turned out 
that interspecies pairing plays the dominant role. To- 
day quantum matter optics offers a means to explore 
condensed-matter and nuclear physics by proxy, such as 
the pairing of fermions in atomic traps and nuclei |5(. 

Here we investigate Raman photoassociation 0,0,00 
of Bose-Fermi mixtures of atoms |l0( , and the subsequent 
prospects for Cooper pairing between different chemical 
species (i.e., atoms and molecules). First, we demon- 
strate that an arbitrary number of stable Fermi molecules 
can be created via fractional stimulated Raman adi- 
abatic passage (STIRAP [L^l ■which is feasible because, 
contrary to bosonic systems 8], collisional interactions 



can be negligible. Density fluctuations in the conden- 
sate leftover from the photoassociation process then re- 
place the vibrating ion lattice of the superconductor |13| , 
and the subsequent phonon exchange can enhance the 
intrafermion attraction. We find that a typical attrac- 
tion is enhanced, but this enhancement is insufficient for 
a transition to atom-molecule Cooper pairs within reach 
of present ultracold technology. 

We model a Bose-Fermi mixture of atoms coupled by 
heteronuclear photoassociation to electronically-excited 
Fermi molecules, which is favored over homonuclear tran- 
sitions for well resolved resonances The excited 
molecules are themselves coupled by a second laser to 
electronically stable molecules. For a degenerate sys- 
tem, the bosonic [fermionic] atoms of mass mo [m+] are 
represented by the field ipo(r,t) [tp + (r,i)], and the ex- 
cited [stable] fermionic molecules of mass m e = mo + m + 
[m_ = m e ] by the field vb e (r,t) [ip~(r,t)], with the boson 
(fermion) field operator obeying commutation (anticom- 
mutation) relations. The Hamiltonian density for said 
non-ideal system is Ti = H.Q + Tir , where 



7~tg 



h 



nv 2 



- [!„ + AocrVo^O 



4>o 



(la) 



\ [(/C + ^V+VA) + n_V» e ^l) + H.C.] . (lb) 



'Presently at School of Pure and Applied Physics, University of 
KwaZulu-Natal, Durban 4041, South Africa. 



The light-matter coupling due to laser 1 (2) is K,+ 
(£!_), and the intermediate (two-photon) laser detun- 
ing, basically the binding energy of the excited (sta- 
ble) molecular state relative to the photodissociation 
threshold, is S (A). Particle trapping is implicit to 
the chemical potential h\x a [a = 0,e,±), and explicit 
traps can be neglected for most practical purposes. 
Low-energy (s-wave) collisions are accounted for by the 
boson-boson (boson-fermion, fermion-fcrmion) interac- 
tion strength Aoo = 2irhaoo/mo (Ao± = / tirhao±/mQ ± , 
A^ = Airfm^ /m* ), with a fflff2 the s-wave scattering 
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length and m* ai(72 the reduced mass. Spontaneous de- 
cay, included as = — T, is generally large enough to 
justify the exclusion of excited-molecule collisions. 

Consider now the process of stimulated Raman adia- 
batic passage from atoms to stable molecules. The key to 
STIRAP is counterintuitive pulse timing, meaning that 
the two lasers are adjusted so that in the beginning, 
when the population is in the initial state (atoms), the 
strongest coupling is between the intermediate and final 
states (excited and stable molecules), while in the end, 
when effectively everything is in the target state, the cou- 
pling between the initial and excited states is strongest. 
As the system passes from atoms to stable molecules, the 
state with the larger population is always weakly cou- 
pled to the electronically-excited molecular state, and 
the subsequently low (ideally zero) population reduces 
(eliminates) radiative losses. 

An overview of our STIRAP model is now presented 
(see also Appendix). When the number of bosons is much 
greater than the number of fermions, Nb 3> Np, the 
frequency scale for atom-molecule conversion is set by 
^+ = \/Pb^+- Although a maximum Nb = 100 is used, 
qualitative scaling to large boson number is cinched by 
assuming a density, ps = 5 x 10 14 cm -3 , consistent with 
Nb = 1.3 x 10 6 Bose-condensed 87 Rb atoms in a trap 
with radial and axial frequencies lu t . = 100 x 2tt Hz and 
u) a = 10 x 2tt Hz. For this density, a ballpark peak value 
for the atom-molecule coupling is S1+ ~ f2o = 27rMHz. 
A typical spontaneous decay rate is T = 10 x 27rMHz. 
The (assumedly negative) Fermi atom-molecule scatter- 
ing length is estimated as |o_| | = aoo = 5.29 nm. The 

number of fermions is restricted to Np = 4 for numerical 
ease, and large-particle-number scaling is again ensured 
by picking a density, pp = 1.1 x 10 12 cm -3 , consistent 
with N F = 5x 10 3 40 K atoms in the same (mass-adjusted) 

trap as the bosons, so that A^ = /V = 5.81 x 27rHz. 

Numerics are further eased by restricting collisions be- 
tween fermions to a bare minimum: ki + k2 = k : } + k4j 
ki = k3, k2 = k4. Also, Nb Np means that collisions 
with condensate atoms can be neglected. 

Of course, complete conversion would leave no atoms 
to form Cooper pairs with molecules, so we pursue 
fractional STIRAP [ll| via the Vitanov et al. pulse- 
forms: n + (i) = fio sinaexp[— (t — t) 2 /T 2 ] and f2_(i) = 
ft exp[- {t + r) 2 /T 2 } + n Q cosaexp[-(t-T) 2 /T 2 }, 
where tana sets the final population fraction. Using 
Fock states and the Hamiltonian given in the Appendix, 
we solve the Schrodinger equation numerically. The 
key results are presented in Fig. ^ For Nf = 4 and 
Nb = 100, the system reproduces the single particle case 
(Nb = Np = 1), i.e., the results are identical to those for 
a three-level atom, as expected for a mostly-undcpletcd 
boson field. Atom-molecule collisions are negligible for 

couplings as large as A_| /fio = 10~ 5 , which is similar 

to a bosonic system, lending confidence to the restricted- 
collision model. Furthermore, Np = Nb = 4 shows 
that many-body effects can limit molecular conversion. 
This many-body effect is similar to the one-color BEC 
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FIG. 1: Creation of an arbitrary mixture of Fermi degen- 
erate atoms and molecules via fractional stimulated Raman 
adiabatic passage. Units of frequency are set by the choice 
f2o = 27rMHz = 1, and the pulse parameters are a = 0.147T, 
T = 5 x 10 3 , and r = 0.7T. For iV s = 100 and N F = 4, 
fractional STIRAP exhibits no visible difference from the 
Nb = Nf = 1 case, while Nb = Nf = 4 illustrates that 
many-body effects can limit the conversion efficiency. 



case [15], and is not attributable to fermion statistics. 

It is not entirely surprising to find a fermionic atoms 
responding as a unit to form fermionic molecules: co- 
operative fermionic behavior w as p redicted for one-color 
formation of Fermi molecules [la ], as well as in atomic 
four-wave mixing |17| and, more relatedly, cooperation 
was observed in magnetoassociation of Fermi atoms into 
Bose molecules 18]. However, what is surprising is that 
STIRAP should work basically as expected even in the 
presence of collisions. This situation arises because, 
when Nb S> Nf, the condensate density is effectively 
fixed, and the associated mean-field contribution simply 
amounts to a static bias that can be absorbed into the 
detuning; with Fermi-Fcrmi collisions blocked, which is 
most likely for small ground-state molecules, only colli- 
sions between the Fermi atoms and molecules can move 
the system off the required two-photon resonance, and 
STIRAP works better. In other words, we get the ad- 
vantage of Bose enhancement of the free-bound coupling 
(f2+ oc J~p~b), without the mean-field shifts. While lim- 
ited computational resources preclude explicit investi- 
gation, these results ought to scale qualitatively with 
increasing particle number, as we have seen for up to 
Np = 20 in one-color production of Fermi molecules. 

Now we are safe to presume the existence of an 
arbitrary admixture of Fermi-degenerate atoms and 
molecules, and thus to consider any subsequent anoma- 
lous pairing. Once the transient STIRAP pulses have 
vanished, the system is described by Ho [Eq. ((Taj) ] with 
A = S = and a — 0,±. For equal-mass fermions, it 
is known that a fermion density fluctuation gives rise to 
an effective chemical potential for the bosons, which ere- 
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ates a bosonic density fluctuation, which in turn leads 
to an effective chemical potential for the fermions. In 
other words, phonons spawned by BEC density fluctua- 
tions are exchanged between the fermions, altering their 
interaction. Just like lattice vibrations that drive the 
attraction between degenerate electrons in superconduc- 
tors, BEC density fluctuations lead to an attractive in- 
teraction that can enhance overall attractions, and thus 
Cooper pair formation [13j . 

Here the effective Fermi-Fermi scattering length is 
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where H = ln(l + x 2 )/x 2 with x = hkp /niQV s and 
v s = (ps^-Aoo/mo) 1 / 2 is the speed of phonons in BEC; 
hkp <C rn,QV s implies H ~ 1. In other words, the 

effective scattering length can be written /a_| = 

1 + Vff — Vfb, where rjpp (t]fb) is the contribution to 
atom-molecule interactions from fermion-fermion (boson- 
fermion) fluctuations. Implicit to expression J2J is the 
perturbative assumption t\fb *C r/pp. The immediate 
contrast with Ref. [13| is that t\pb < is allowed. For 



a weakly attractive system {kp\a aia2 \, ps^axa^ 
the critical temperature for Cooper pairing is 
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T c = 0.61 T F exp 



tt/4 



(3) 



where Tp 
with Rm : 



h{p+ + /i_)i?M /ks is the Fermi temperature 



was taken as the same for both species, so that = 
(m±/m* )p.± and Tp = Tp -J m + /m_ . 

Continuing to focus on 87 Rb- 40 K, N B = 1.3 x 10 6 
BEC atoms in a trap with u r = 100 x 27rHz and 
u> a = 10 x2ir Hz yield a boson density ps — 5xl0 14 cm~ 3 . 
A modest number of fermions, Np = 5 x 10 3 , means that 
the loss of condensate atoms in molecule formation can be 
neglected, the condensate will absorb any heat created by 
pairing- induced holes in the atomic Fermi sea |l9j . and 
collapse instabilities (Modugno et al. 0) are avoided. 
Presuming that fractional STIRAP converts roughly 18% 
of the initial Fermi atoms into molecules (see Fig. 1), and 
that the fermions share the same (mass-adjusted) trap, 
then p± = 1.1 x 10 12 cm -3 , and the requirement of equal 
Fermi wavenumbers for the atoms and molecules is met. 
For the given parameters, the size of the BEC is roughly 
an order of magnitude larger than the Fermi clouds, so 
that overlap should not be an issue. 

Figure [21 summarizes our investigations. Under the 

above trap conditions, and for |a_| | = ago = 5.29 nm, 

we find the weakness parameter kp\a^ | ~ 0.05 and 

an unreasonably low critical temperature. Nevertheless, 
if the trap is modified to be anharmonic, then tighter 
confinement ultimately means a diluteness parameter on 

the order of kp\a^ | = 0.2 |2jj, and the situation is 

improved but still out of reach of current technology 
(T ~ 0.057> [3). The only chance appears to be for a 
tight trap and a large scattering length, kp\a^ | = 0.5; 
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FIG. 2: Critical temperature for the superfluid transition to 
anomalous atom-molecule pairs as a function of fermion-boson 

fluctuation strength. Calculations are for weak (fcir|a-i | = 

0.05), marginally weak (0.2), and marginally strong (0.5) in- 
teractions. The solid (dashed) curve gives the critical tem- 
perature in units of the effective (atomic) Fermi temperature. 



however, the theory is at best marginally applicable in 
this regime and, besides, some other process (e.g., three- 
body recombination) would probably win out before su- 
perfluidity could set-in to such a system. 

In conclusion, collisions can arguably be made neg- 
ligible in photoassociation of Bose-Fermi mixtures of 
atoms, and fractional stimulated Raman adiabatic pas- 
sage scheme is therefore feasible for creating Fermi de- 
generate molecules. On the other hand, while possible in 
principle, accessible atom-molecule superfluidity means 
that next-generation technology must shed another or- 
der of magnitude in temperature from the first gener- 
ation 21], or that a system with an attractive and a 
strong-but-not-too-strong interaction will be found. Our 
opinion is that it will be worthwhile to look for other 
ways of driving atom-molecule Cooper pairs, such as a 
photoassociation or Feshbach resonance. 

Appendix.-To model STIRAP from Eq. make a 
time-dependent unitary transformation U(t) — Ii a U a (t) : 
where f7 0(+) (t) = exp[-ifi +(Q) ^ +{Q) ^ +(Q) t}, U e (t) = 

exp[i(^ + + n )iplip B t], U-(t) = exp[-i(^ + + /j, )ip'_ij}^t, 
whereby fj, e — > p:' e = /i e — (/i_ + po) and — > p!_ 
ji — (jj,_ + fig); in turn, absorb A*e(-) m t° t ne intermedi- 
ate (two-photon) detuning. Assume that atom-molecule 
conversion occurs on a timescale much faster than the 
motion of the atoms in the trap, and thereby neglect 
the kinetic energies and any explicit trapping potentials. 
Take the Fermi energy to lie within the Wigner threshold 
regime, so that the coupling /C + is the same for all modes. 
Focus on the regime Nb ^> Np, so that the condensate 
is practically undepleted and ipo ~ \fp~B can be absorbed 
into the atom-molecule coupling strength; likewise, the 
BEC mean- field shifts, Xoaw^o ~ PB^Oa, are a constant 
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that can be absorbed into the detuning, and therefore 
of no concern. Input two-photon and intermediate reso- 
nance by setting A = and S = — hiT. Finally, switch 
to momentum space, so that 

T = ~5 ir XXe ake + A +- a L + a k 2 - ak 3- a k4 + 

k {kl} 

[(^+ a L a k+ + ^-aLak-) + H.c. , 

k 

where VL+ = ^/pbIC+ and = A_| jV . Recover the 



few-boson results by assuming that Bose stimulation fa- 
vors condensate modes, and making the substitutions 
fl + — » n+fy/Ns and aj ce ai c + — > a^a^+aQ. 
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